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Under the assumption that a wing, body, or wing-body combination 
is slender or flying at near sonic velocity, expressions are given which 
permit the calculation of pressure in the immediate vicinity of the con- 
figuration. The disturbance field, in both subsonic and supersonic 
flight, is shown to consist of two-dimensional disturbance fields 
extending laterally and a longitudinal field that depend on the stream- 
wise growth of cross-sectional area. A discussion is al given of 
couplings, between lifting and thickness effects, that nec 3sarily arise 
as a result of the quadratic dependence of pressure on the induced 
velocity components. 


INTRODUCTION 


This paper is concerned with the prediction of pressure distribu- 
tion on or in the immediate vicinity of a wing, body, or wing-body com- 
bination under conditions in which the geometric configuration is 
slender in the flight direction or is flying at near sonic velocity. 

The material to be presented is thus associated with the rather exten- 
sive group of results that belong to what is often referred to as 
slender-wing theory. The basic assumptions and methods can be found in 
publications by Munk, R. T. Jones, and Ward (refs. 1,. 2, and 3) 1 and. a 
discussion of the applicability of the methods to the prediction of 
loading on slender wings at sonic flight speeds has been given in refer- 
ence 5. In reference 2, attention was directed toward the calculation 
of load distributions over wings in subsonic and supersonic flight and 
reference 3 was devoted essentially to the consideration of supersonic 
flight velocities. It is therefore of interest to investigate further 
the effects attributable to thickness on wings and wing-body combina- 
tions at both subsonic and supersonic flight speeds. Such investigations 
lead to valid approximations of interference effects and also indicate 
the way in which thickness and lifting effects can produce couplings in 
the calculations of pressures induced in the flow field. 

Reference should also be made to the recent extensions of slender-wing 
theory by Adams and Sears (ref. 4) . 
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ANALYSIS 


It is proposed to take the basic solutions of the linearized 
partial differential equations governing three-dimensional compressible 
flow and to obtain a simplification of the expressions by restricting 
attention to the induced field in the immediate vicinity of slender 
airplanes or missiles. These simplified expressions contain solutions 
used previously to study the forces and moments on lifting wings and 
bodies. In addition, however, they can be used to evaluate the first- 
order thickness effects on the pressure in the vicinity of the wing 
and body. 


Consider, first, the construction of a weakly disturbed flow field. 
Let a uniform stream flow in the direction of the positive x axis of 
a Cartesian coordinate system, as in sketch (A) . Immerse in the stream. 


A 



•vdiich has a velocity U 0 and a 
Mach- number Mq, a slender wing- 
body shape the surface of which is 
inclined at a small angle to the 
free-stream direction. This angle 
of inclination must be small enough 
so that nearly everywhere in the 
fluid the magnitude of the pertur- 
bation velocity vector divided by 
the speed of the free stream is 
much less than one; that is, 

£EI3l±2. <<x ( ia) 


Moreover, large supersonic Mach 
numbers are to be avoided and as a 
measure pf this condition the 
inequality 

Mo ^ u2 + v 2 + J g <<i 

Uo 

is imposed. 


Consider, next, the linearized partial differential equation gov- 
erning weakly disturbed isentropic fluid flow. In terms of the pertur- 
bation velocity potential <p(x;y.z), the lowest order approximation 
consistent with inequalities (la; and (lb) is 


(1-Mq 2 ) 9^+ <Pyy +<P 2Z = 0 


( 2 ) 
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where the subscripts denote partial differentiation with respect to the 
indicated variable. 

Consider, finally, the expression for the pressure coefficient 
that is again consistent to the lowest order with inequalities (la) 
and (lb) . By expanding the pressure -velocity relation for steady 
isentropic flow and neglecting’ higher-order terms, one finds 

2u (1-Mo 2 )^ + v2 + w 2 
T 0 U? ~~~ 

where p and p are pressure and density, respectively, and the sub- 
script o refers to conditions in the free stream. It follows from 
inequalities (la) and (lb) that pressure coefficient can be expressed 
in the form 


Cp = 


P-P, 


ibo u o £ 


ft! — 


Cp 


ss 


2u 

T 0 


V 2 + w 2 



(3) 


Equation ( 3 ) is the simplest general expression for pressure coefficient 
that is still entirely consistent with the assumptions basic to the 
development of equation (2). 

Special solutions applying to problems of the class indicated can 
be obtained by appropriate simplification of general solutions to 
equation (2) . Such a procedure will be discussed in the next section. 
The pressure coefficient is then determined by substituting these 
results into equation ( 3) • The simplifications that can be made in 
evaluating the pressure on the surface of the airplane will also be 
discussed. 


The Reduced Solutions 


Subsonic . - As it applies to subsonic flow, equation (2) can be 
written in its normalized form as 

cp + cp + cp = 0 (^) 

Rxx ^yy T zz * 

The analysis of equation (*J-) can be interpreted as applying to the condi- 
tion Mq = 0 but one can extend the solutions throughout the subsonic 
Mach number range by applying the Prandtl-Glauert rule . It is important 
to stipulate, however, that the term "slender” will, unless otherwise 
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indicated, have a dual interpretation - describing an aerodynamic con- 
figuration that is either much longer than it is wide or is flying 
close to the speed of sound. 


A well-known solution to equation (4), resulting from an applica- 
tion of Green's theorem, is given by the expression 


<p(x,y,z) 



dS x 

y (x-xx) 2 + r 2 


(5) 


where dS x is the eleme nt of surface area on the airplane or its 

vortex wake, r equals J (y - y^) 2 + (z - z-j^) 2 , and bfdn' is the 
derivative normal to the surface S. When this solution is applied to 
boundary -value problems for slender configurations it can be simplified 
considerably. One method for bringing this about is suggested by 
studying the variation of [r/(x - x x )] 2 over the area S. Consistent 
with the assumptions made, one has [r/(x - x x )] 2 < < 1 over almost all 
-of the airplane surface and vortex wake provided the point x,y, z is 
on or in the vicinity of these surfaces. This implies the approximation 



Xi ) 2 + r 2 ~ |x 



( 6 ) 


The singularity at x = x x , which thus appears in the integrand of 
equation (5) for the limiting case r = 0, produces a divergent integral. 
An indication of the manner in which this difficulty can be avoided is 
obtained through consideration of the single integral 



g(x 1 )dx 1 
.Ax-x,) 2 + r 2 


If g(x z ) is differentiable within the region of integration, it 
is easy to show that this expression can be written in the form 


a 1 * , x , ^-1 x " x i „ a r* r \ x-xi , Ix-xxi+ytx-x,)^ 

-J g(x,) slBh g WKi ! " r — 


and the approximation for small values of r then becomes 

■>b x-x. 2|x-x x j 


^ po x-x, 

~ / s ^ 1 ) 1 V n 

c)x J 1 |x-x x 


dx. 
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Equation (5) is to be -written subsequently with a logarithmic 
kernel in the integrand and a derivative operator outside the integrals. 
An estimate of the order of accuracy involved in rising the approximate 
form of the integrand follows from the evaluation of the two expressions 

i 0 = r ava*,. 

J Bo J 

and 



Sketch (Bl) shows the variation of the ratio E = (l 0 -I 1 )/lo f° r a 
rectangle wing of area S Q and aspect ratio A = l/5 as the point x,y 
covers the portion of the xy plane on and within one semispan of the 
wing. 


The maximum value of E occurs when 
the point x,y lies along the trailing 
edge (or, by symmetry, along the leading 
edge) and it is significant to notice that 
the value of E decreases as the point 
moves from this location in either x 
direction and increases as it moves from 
the x axis in either y direction. 

This illustrates the necessity of restric- 
ting the approximate solutions to slender 
configurations and, further, to portions 
of the flow field in the vicinity of the • 
longitudinal axis .of these airplane 
shapes. Sketch (B2) shows how the maxi- 
mum value of E within one semi span of 
the wing decreases with decreasing 
aspect ratio. 

Under the restrictions that have 
been imposed, it is justifiable to 
introduce simplifications in the form of 
the derivative d/chi' and the differ- 
ential area dS x appearing in equa- 
tion (5) • The operator d/dn' can be 
expressed as 

nj_d n^ 
dx dy dz 
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•where e-ijQ -2 n.nfl hq are the direction cosines between a normal to the 
surface S and the x,y, and z axes, respectively; the differential 
area dS x can he expressed as 


ds-j^dx-L 

\/l - n x 2 

where ds x is a differential length along the surface in a yz plane. 
If the airplane is slender, n x is small, and can he neglected, rela- 
tive to either unity or ,/n 2 2 + n 3 2 . 

Combining the two simplifications discussed above, one can approxi- 
mate equation (5) by the expression 




JL_ _c*_ 
tit dx 



x-x x 2|x-x x 

; t ln ~ 

X-X, r 


(T) 


where S fdn represents n 2 S/Sy + n^fbz f the normal derivative to a 
section in a yz plane, s is the curve bounding this section, and Z 
is the total length of the airplane. Equation (7) reduces to 


q>(x,y,z) 




i a r tel -as- 

4lf Bx J _ x-x x 



f s) lir 

ds x 



Since 



s 


(8) 


where S(x) is the cross-sectional area of the specified shape in a yz 
plane, this expression becomes 


cp(x,y,z) 



jr f s ' ( x i)f— — in 2 (xi-xjax 

tir ox J |x-x x 


It is apparent that the first integral on the right-hand side of 
this equation is a solution to Laplace’s equation in two (the y and z) 
dimensions. This portion of the solution will be represented by the 
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symbol <p 2 (x;y,z). The x dimension, -which does not appear explicitly, 
enters as a parameter when the solution is adapted to particular 
boundary values- Hence, b y means of the usual Prandtl-Glauert trans- 
formation x — >x, y — >y a/i-Mq 2 , z — > z a/i-Mq 2 , together with equa- 
tion (8) and the relation 

h{ s ' (Xl) l^F^ - ' 2S ' (x) 

the final form of the solution is 


qp(x,y,z) = cp 2 (x;y,z) s ' ( x i) 


x-x. 


In 


2 


x-x. 


l X " X l| a/Am^ 


dXi 


(9) 


The physical significance of the last expression can be interpreted 
easily. The three-dimensional perturbation velocity field induced by 
wing-body shapes that are slender or flying at near sonic speeds is 
approximated in the vicinity of such shapes by a field that satisfies 
the two-dimensional Laplace equation and the boundary conditions in 
transverse planes plus a longitudinal field which depends on the stream- 
wise rate of change of cross-sectional area and is independent of y 
and z. 

Supersonic. - In the case of supersonic flow, the normalized form of 
equation (2) becomes 


^xx “ 'Pyy “ ‘Pzz = 0 


(10) 


An analysis based on equation (10) applies specifically to the condition 
Mo = a '/"2 but these results can be extended throughout the supersonic 
Mach number range by applying the Prandtl-Glauert rule. Vol terra’s 
solution to equation (10), (see, e.g., ref. 6, p. 190) which is analo- 
gous to the subsonic form given in equation (5), is expressible as 


<P(x,y,z) = - 


A. A 

2jt 8x 



T 



Zn 


X-X! 


+ 


J (x-Xj) = 


dSi (11) 


where, as in the subsonic case, dS^ is an element of surface area on 
the airplane or its vortex sheet and r equals «/(y - y x ) 2 + (z - 
In distinction to the subsonic solutions, the area of integration is 
now the portion of the airplane and its vortex wake within the forecone 
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from the point x,y,z, and bfbv is the derivative along the conormal 2 
rather than the normal. 

If the application of equation (ll) is limited to slender configu- 
rations, the approximation, similar to expression (6) for the subsonic 
case, 

sa jx-xxl 


is implied. Furthermore, if the conormal and differential area are 
expressed in terms of the directi on cosines and n x is again neglected 
relative to unity or J n 2 2 + n 3 2 , equation (ll) reduces to the form 


<P = 


2rt dx 




Zn rdsx - 


_1 _ 8 _ 
2ir 8x 




89 

8n 


dS! 

( 12 ) 


This differs from the corresponding expression, for subsonic flow only by 
a factor of two and by the extent of the x x integration. The latter 
is carried only to x in the supersonic case since the original inte- 
gration area r included only points in the forecone from x,y,z. It 
is obvious, however, that these two differences are compensating in the 
first integral term so the final expression for the perturbation poten- 
tial, appropriately modified by the Prandtl-Glauert rule, becomes 


<P(x,y,z) 


<P 2 (x;y,z) 



2(x-xi). 
/ Mo 2 -1 


( 13 ) 


The Reference Coordinate Systems 


Equation (2) was developed specifically for the case in which the 
undisturbed stream at infinity is parallel to the x axis. A coordinate 
system so orientated is usually referred to as the wind axes. (See 
sketch (c).) When the configuration is tilted with respect to the 


2 The conormal is the vector that results from changing the sign of 
the x component of the normal. 
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free-stream vector, however, it is often easier to study the boundary- 
value problem with the x axis placed along the center line of the 
fuselage . Such a coordinate system is usually referred to as the body 
axes. 


Obviously the wind and body axes differ significantly only by rota- 
tions about the y .and z axes. When Mq is zero, equation (2) is 
invariant to such a rotation-, but for values of Mo greater than zero 
this is no longer true. However, -when Mq is greater than zero, equa- 
tion (2) represents the governing differential equation only to a 
certain order, and, if the magnitude of the rotation is similar to that 
of the parameters by which the equation is ordered, it is, in this sense, 
still invariant to rotations about all three axes for both subsonic and 
supersonic Mach numbers. Thus equation (2) is to the lowest order the 
governing partial differential equation for both wind and body axes, 
provided the airplane is slender and the angles of attack and sideslip 
are small. 

Although the partial differential equation is invariant with respect 
to a small rotation of the coordinate system, the boundary conditions and 
expression for the pressure coefficient in terms of the perturbation 
velocities are not. The following will contain a 'discussion of the 
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boundary conditions and tbe pressure and loading coefficients with, 
reference to a body axes system. 


Tbe Boundary Conditions 


The boundary conditions require that the gradient of the total 
velocity potential evaluated infinitely far from the aircraft be con- 
sistent with a uniform free stream there (the direction of "which 
depends on the orientation of the coordinate system) and "when evaluated 
normal to and on the surface of the airplane itself be zero. Let 
$(x,y,z) denote total velocity potential, cp(x,y,z) perturbation velocity 
potential, and refer the analysis to body axes in a free stream. If 
the orientation of the free-stream velocity vector to the system of axes 
is fixed by the angles a and 7 as shown in sketch (C), one can write 


<l>(x,y,z) = U 0 (x cos a cos 7 +y sin 7 +z cos 7 sin a) + cp(x,y,z) 


such that on the aircraft surface 

Uo(na. cos a cos 7+n 2 sin 7+113 sin a cos 7) + rLx<T> x + ^z^y + = 0 


n i> n 2 , an< l H 3 again being the direction cosines of a normal to the 
airplane surface with respect to the x,y, and z axes, respectively. 
By the assumptions basic to the. present theory, the latter equation 
reduces to 


U Q (n 1 +n a 7 +n 3 a) + — <P £ (x;y,z) = 0 


(Ik) 


where, as before, n is the normal to the curve bounding a cross section 
in the yz plane. 

Equation (l4) , which applies to arbitrary slender shapes, can be 
simplified for many specific problems. Consider now three types of 
configurations that lead to such simplifications: first, a surface, 

such as a wing, which deviates only slightly from a plane; second, a 
surface which forms a body of revolution; and, third, a surface which 
is a combination of the above two. 


Planar problems. - Let h(x,y) be the distance a surface deviates 
from the z = 0 plane, and b represent the wing span. (See 
sketch (Dl) . ) Assume that a/(db/dx)<<l holds; then furthermore. 
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if the inequality (dh/Bx)/(db/dx) < < 1 is satisfied, it is consistent 
with the previous approximations to ' 


neglect the y component of the 
normal along the wing surface and 
to project the yelocity vector 
represented by the resulting verti- 
cal derivative to the upper or 
lower surface of the z = 0 plane. 
In this way equation (1^) becomes 


U 0 (n 1 +n^x)+ n 3 



= 0 


and, since nx/n 3 = -cSh/dx, the 
boundary conditions for planar 
problems 3 are expressed by the 
equation 



Body of revolution problems . - 
Let B be the radius of a body of 
revolution. (See sketch (D2).) If 
9 is measured from the z axis in. 
the yz plane and 7 is set equal 
to zero, the relations 

n l 




together with equation (l4) , give 
for the normal derivative on the 
surface of the body 



U 0 — — Una cos 9 (16) 

dx 



Sketches (Dl), (D2), and (D3) 


Certain planar problems, such as the cruciform wing, require more than 
one plane but the concepts are essentially the same as those presented 
here . 
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which is the simplified expression of the houndary conditions pertaining 
to bodies of revolution. 

Interference problems. - Consider, finally, surfaces which are a 
combination of the above two as, for example, the one shown in sketch (D3) • 
The rather obvious extension of the above concepts is to apply equa- 
tion (15) over the winged portion of the configuration and equation (l6) 
over the body. It is then necessary, however, to consider the relative 
magnitudes of the terms' 8h/8x, ds/dx, dt/dx, and dR/dx, since they 
appear in the solutions in various combinations. If the winged portion 
is to be treated as a planar problem, the magnitude of cJh/8x must be 
small enough to be neglected in comparison to the leading- and trailing- 
edge slopes, ds/dx and dt/dx. But this does not imply that 8h/8x can 
be neglected in comparison to dR/dx or that dR/dx can be neglected 
in comparison to either ds/dx or dt/dx. The latter approximations will 
not, in general, be made. 


The Pressure Coefficient 


The expression for the pressure coefficient given by equation (3) 
is written in terms of velocity components that are referred to the 
wind axes. Its re-expression in terms of velocities referred to the 
body axes is readily determined. For the orientation shown in 
sketch (C2) , the equation becomes 


Op <?> x + Wy + <%) OP/ + *■?) (17) 

Equation (17) can be -used, in general, to evaluate the pressure in 
a perturbation velocity field that is referred to the body axes. If 
the interest is limited to the pressure on the surface of the aircraft, 
however, certain simplifications can be made. For example, consider 
the configuration illustrated in sketch (D3) consisting of a swept-back 
wing mounted on a body of revolution. For simplicity, let 7=0. 
Applying the boundary conditions given by equations (15) and (1 6), one 
can show that on the surface of the wing 


°P =- 


- U 0 


+ 




and on the surface of the body 


C p =- 


29. 


L u 0 


^ + _i_ (U D ct sin 0- i cp e )' 


Ur 


R 


+ a 2 


B 



(18a) 


(l8b) 
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These solutions can "be simplified further "by considering the 
detailed nature of the perturbation velocity field induced by shapes 
such as that shown in sketch (D3) • Thus, the results given by equa- 
tions (9) and (13) can be expressed in the form 

qp(x,y,z) = <p 2 (x;y,z) + A(x) (19) 


where the expression for A(x) depends on whether the speed is subsonic 
or supersonic. Further, for the particular configurations being con- 
sidered, the expression for <p 2 (x;y,z) can be written in the general form 


cp 2 (x;y,z) = a<P a (t,s,Rjy,z) + <|^(t,s,Kjy,z) + ~ <P c (t,s,R;y,z) (20) 


dx 


since the dependency on x can enter only through the boundary condi- 
tions which, in turn, 'are specified by the body radius R(x), the wing 
thickness, h(x,y), and the lateral distances from the center line to 
the trailing edge and leading edge, t(x) and s(x), respectively. The 
term. acp a will be referred to as the potential due to angle of attack, 
since it vanishes when the angle of attack vanishes and increases 
linearly with increasing a; the term (Bh/8x)cp^ + (dR/dx)q> will be 
referred to as the potential due to thickness, since it exists idaen the 
angle of attack is zero, does not change with angle -of -attack change 
and vanishes when the thicknesses of the wing and body do not vary 
with x. 


By breaking <p 2 down into its component parts as in equation (20), 
it has been ordered in that the magnitudes of the terms on the right- 
hand side of equation (20) are controlled by the coefficients of the cp’s, 
and the derivatives of cp a , cp^, and cp c with respect to s,t,R,y, and z 
can all be considered equal. Since a and Bh/Bx are negligible relative 
to dt/dx, ds/dx, and dR/dx (as was- pointed out in the discussiop of the 
boundary conditions for interference problems), equations (l8a) and (l8b) 
can be written:, 

on the surface of the wing 

(21a) 

By By ox By By dx \ By/ J 9 Z _ 0 


f 2_ + _1_ dR 

lu 0 Bx u D 2 <3x 
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and on the surface of the body 

2 B<P 


= 2_ Va E\ 2 + __i M 2a ^a [2 Bh ^c , 

lu 0 Bx \.dx/ Uo^ 2 Se Be a* Be Be 


® fdV ^‘ - 2o0oE sin 0 ^ 


dx V Be 


Be J 


(21h) 


If the hody is a cylinder so that its radius does not vary -with x, 
the pressure coefficient reduces to 


Cp - 


_2_ Bq> 

L TJ 0 Bx J 


( 22 ) 


B 


Loading Coefficient 

By definition the loading coefficient is 


£P 


— = (Cp) z - (c p ) u 


(23) 


■where the subscripts u and l refer to the upper and lower surfaces of 
the airplane, respectively. It is immediately apparent from, an inspec- 
tion of equations (23) and ( 19 ) that the loading is not affected by A(x). 
Hence, the lift, pitching moment, rolling moment, and induced drag can 
all be expressed entirely in terms of cp 2 (x;y,z). 

Consider again the type of airplane shapes represented in 
sketch (D3) . The velocity potential cp £ for such a class of conf igura- 
tions has been expressed in equation (20) as the sum of three potentials: 
one due to angle of attack, one due to the thickness of the wing, and 
one due to the thickness of the body. It is now useful to remark that cp a 
has odd symmetry with reference to the z = 0 plane and cp^ and cp c 
have even symmetry. Placing equations (21a) and (21b) into equation (23) 
and using these properties, one finds 

(W\ _ 2a dq> a 1 dR ^q> a Bq> c \ 

' 3- 'wing |_ u 0 dx V By By / 

z=o 


(24a) 
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and 


*4 


ody 


2a 


+ 1 

®a! 


d<Pc\ " 

U 0 

L 3x 

UqR 2 

dx 

v Be 

Be/ - 


(2te) 


B 


where A indicates the difference between a quantity on vertically 
opposed points of the tipper and lower surface of the airplane. 


It is apparent from the last two equations that, in general, the 
angle -of -attack and thickness solutions have a coupling effect on the 
loading coefficient and therefore their contribution to the load distri- 
bution cannot be treated separately. It is also important to notice the 
two special cases in which the coupling effects vanish; namely, a body 
of revolution without wings, and an airplane with a cylindrical body 
between the foremost and rearmost extent of the wing. In the former case 
the term dcp c /ci 6 is zero and in the latter dR/dx is zero. In both 
these cases the equation for the loading coefficient is 


SE . **3 . J. A 

i- Bx U D Bx 

The Total Lift 


( 25 ) 


Total lift can be obtained, of course, by integrating the loading 
coefficient over the aircraft surface. A much simpler way of finding 
the lift, however, can be derived from a momentum balance. Thus, by y 
momentum considerations it is possible to show that the vectorial force F 
on a body inside a control surface S is given by the surface integral 


-> p p ->■ npf->-^\->-> 

F = -J J (p-Po)dS - J J p\V - Y 0 J |_V- dS_ 

S S'. 

where vector notation is used, the o subscript indicates free-stream 

conditions, p and p are the local static pressure and density, and Y 
is the local velocity vector. Let the surface S be a cylinder of 
infinite radius and two yz planes closing the cylinder and located 
infinitely far ahead of and behind the airplane. Then the lift force 
is given to the lowest order by 


L 



[pwCUo+u)]^ dy dz 
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which reduces to 


00 00 

l‘-t>cPof f (v) x „ 

-00 -CO 


dy dz 


This can he simplified since w = Scp/Sz and (A9) x=00 is the same as 
the jump in the potential evaluated at the airplane trailing edge. 

Thus the expression for lift becomes 

(26) 


Equation (26) applies to all slender shapes. In special cases 
represented by sketch (D3) , thickness effects always have even s ymme try 
with respect to the z = 0 plane, and it follows that the total lift 
and the vortex distribution in the wake of such configurations are 
affected only by the part of the potential due to angle of attack, even 
though the detailed load distribution depends upon both thickness and 
angle -of -attack solutions. 


L = P n U, 


°f 


span 


(A<P) T _ E _ ay 


EXAMPLES 

Pressure on a Triangular Wing 
With Elliptic Cross Section 


It is of interest to calculate, by equation (13) , the pressure on 
nonlifting wings of triangular plan form and elliptic cross section 
flying at supersonic speeds, since examples of this type have been 
solved without restriction to slender-wing theory. It is proposed, 
therefore, to study two cases given first by Squire (ref. 7) and then 
to compare the analytical results. 

Let the wing be placed at zero angle of attack in a supersonic free 
stream of Mach number Mq. Consider first the thickness distribution 
for which the ordinate of the upper surface Is 

b(x,y) = — V'n^x 2 - y 3 

2mc 0 


( 27 ) 
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where c D is the root chord, t is. wing thickness at x = Co, and m 
is the tangent of the semi apex angle of the plan form. Since the flow 
is supersonic, it is unnecessary to consider closure. 


Since attention is confined to symmetric nonlifting wings, the 
boundary values are planar and are expressed by equation (15) for a = 0 
Further, the solution is given in terms of these boundary values by 
equation (12) since (5<p/chi) z-0 becom es U 0 (dh/c)x) and (A<P) Z=0 - is zero 
by symmetry. Hence, if 3 = J Mq 2 -1, 


~ q>(x,y,0) 

u o 


mx 


1/ |*»lr 


■7! I dy : 


-mx 


_ i__b_ r 

1 2rt bx J 


S' (xi) Zn 


g(*-*i) 

3 


dx n 


It follows from equation (27) that the elliptic section in the 
plane x = x x has major and minor semiaxes equal to mx and tx/2e 0 , 
respectively. The cross-sectional area and the surface slope are, 
therefore. 


S(x) = ittmx 1 2 /2c 0 , bh. fbx *» tmx/2c 0 J m?x 2 -y s 


Since 


mx 



tmxZn | y-y^ | dy z 
2c ov /m 2 x 2 -y 1 a 


mtx mx . 

2^ 111 IT' 


|y| < mx 


mtx , 1 7 + <y/ 7 2 -m 2 x 2 , 

_ — Zn 1 \y >mx 

2c n p 


the expression for perturbation potential becomes 


~ <P(x,7,0) 
Uo 


fSsTL.la 

I 2Co V 

tmx , 
t — 1 - Zn 

[ 2c o L 



|7|<mx 


tac 

3 (|y| + V y2-m2 x 2) 



>mx 
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From equation (22), pressure coefficient in the plane of the 'wing 
is, 



|y|< 


mx 


kx 


P( |y|W y^-m 2 ^ 2 ) 


- 1 


s|y|> 


mx 


(28) 


The pressure distribution on the wing is uniform and off the wing has a 
square-root singularity at the leading edges. Analysis not limited by 
the assumption of slenderness yields for pressure coefficient on the 
wing 


_ K-E 

P " c o -s/l-p^ 2 

where K and E are complete elliptic integrals with modulus A/l-p^m 2 . 
Since for values of the modulus near one the asymptotic relations 

ll- 

K«I n — , E * 1 

pm 

apply, the pressure coefficient in slender-wing theory is seen to be a 
first-order approximation. 

Squire has also considered the wing with ordinates given by 


h(x,y) = - ■ t - '/m 2 x 2 -y 2 
2c 0 m 


( 29 ) 


The lateral section is again elliptic, with semi major and semi minor axes 
equal to mx hnd tx 2 /2c 0 2 . Cross-sectional area and surface slope are, 
respectively. 


jtmtx 3 



8h t 2m 2 x 2 -y 2 

SC - ■ 

8x 2mc 0 2 V m %^-y 2 


S(x) = 


i 
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By direct integration it can "be shown that 

(30) 

Analysis not limited to slender wings yields the expression 


Op _ — [(3 - 3 2 m 2 )K-(4-2g a in s )E ] 

c^fl-fsV 8 ) 

and again the results are in agreement if higher-order terms in pm 
are neglected. 

A study of generalized conical flow fields in linearized supersonic 
theory reveals that the linear pressure distribution in the above problem 
can also be obtained on a wing with thickness specified by the relation 

fc(x,y) = ^jr cosh" 1 ~ 

y 

where k is a constant that can be related to the ma-TlnTum thickness 
ratio (attained along the line mx/y = I.31) of the wing. The cross- 
sectional area and surface slope are, respectively, 

sW = iSHf, . *** 

3 8x m >J m 2 x 2 -y 2 

and pressure coefficient on the wing is 



The latter expression agrees to the -first order in pm with the general 
linearized solution for such a wing presented in reference 8. Slender- 
wing theory thus retains the property of the more general linear theory 
in that a given pressure distribution does not necessarily yield a 
unique thickness distribution. 


Supersonic Drag of Wings at Zero Incidence 


The general expression for the supersonic drag of a slender aero- 
dynamic shape has been derived by Ward (ref. 3) through the use of 
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momentum methods. It is also possible to obtain these results by direct 
integration of the product of pressure and surface slope over the speci- 
fied surface ; the analysis, however, requires rather careful attention 
to orders of integration, when planar problems axe involved. Consider, 
for example, the drag of a wing at zero incidence and with a specified 
thickness distribution z = ±h(x,y). The drag of the wing is express- 
ible in the form 


D — Dg + 2q 



dx dy 


(31) 


where the first term includes possible contributions to the drag that 
result from a finite leading -edge radius of curvature. From reference 9 , 
this drag per unit of span is, in slender wing theory. 


dD e 

<3y 





(32) 


where r n is the radius of curvature normal to the wing leading edge 
and s is the local semispan. If the ordinate of the wing, in the 
vicinity of the leading edge, is 


z u =f(s,y) A /s - y 


equation (32) becomes 


dy 


*q 


2 



From equation (13) the potential of the wing, evaluated in the 
plane of the wing is 


. U o P dzu(x,y ) U Q ,, , 

q)(x,y,0) = — / — in y-yj ay x - ~ S' (x) 

J D £x 2rt 



Uo 

2rt 


f S^Xj) In (x-Xi) dx-L 


(33) 


and, since pressure coefficient in the planar case is directly propor- 
tional to the streamwise gradient of cp, the contribution of each of the 
terms on the right-hand side of equation (33) can be calculated separately 
in equation (31)- The second and third terms offer no difficulty but 
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simplification of tie expression resulting from tie first term neces- 
sitates an inversion of order of integration and, if tie leading edge 
ias a finite radius of curvature, suci an inversion cannot le carried 
out in tie conventional manner. However, a metiod by means of wiici 
suci an inversion can be carried out is presented in reference 10. 
Thus, set 


Ii 



dzufoy) d g Zu(*jyx) 
dx Sx 2 


In 



and 


if . r\ Zzu(x,y) ^ 2 zu (x,y x ) , , 

1 2 = r dy / <ay - — IrL |y-y \ 

where J- refers to a finite -part integral 4 and tie notation J' dy J* dy x 


signifies that tie y 1 integration must be performed first. Then if 


Zu(x,y) = f (s,y) Js - y 


it can be shown that 

11 -^"T (s) 12 M 

Detailed analysis reveals that tie residual term (i.e., tie value of 
II " I 2 ) yields a drag component that is equal in magnitude but opposite 
in sign to D e . 


Tie final expression for tie drag of tie wing is then 


D 

q. 


-»(l) 


(2) ^i u (x,y)] ^ ,r S(Z )a hu(x,y 1 J 


dx 


x=l 


f 


-B(l) 


5x 


x=2 


Injy-ylcLy, + — & ■ In | + 

A x 2it p 


is '(!) J T I S~(x 1 ) in (I-x x ) ax r Xj l S'U)Sxf S*^) Inhc-zja*,. 

0 0 CM 


where 2 is tie over-all length of tie wing. 

^Tor a definition of the finite -part integration technique as used here, 
see reference 10’. 
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As a particular example, consider the wing -like surface of trian- 
gular plan form (ref. 7) which results from a combination of the 
surfaces specified in equations (27) and ( 29 ) and has ordinates given 
by the expression 


ku(x>y) = (e 0 — x)V^W (35) 

mc 0 2 

This wing has rounded leading edges and a finite trailing-edge angle, 
and from equation ( 3 ^) its drag coefficient based on wing area is found 
to he 

Cd = + (36> 

It is apparent from equation (3^0 that wing drag varies with Mach number 
so long as the streamwise gradient of area is finite at the rear of the 
wing; conversely, there is no dependence on Mach number when the gradient 
of area vanishes there. For example, a wing with an elliptic plan form 
and biconvex sections satisfies the latter condition, and its drag coef- 
ficient based on wing area is 

C D = 2 1 (t7 (37) 

where t is total maximum thickness, a is the semiaxis of the elliptic 
plan form in the stream direction, and b is the semiaxis measured 
normal to the stream direction. 

A comparison between the values of Cp given by slender airplane 
theory for the Squire wing (eq.. ( 36 )) and the elliptic lens (eq.. (37)) 
and the exact thin -airfoil -theory values 5 for the same wings is shown 
in sketch (e) . 


5 The exact values were given, respectively, by Squire in reference 7 
and by R. T. Jones in a paper entitled "Theoretical Determination of 
the Minimum prag of Airfoils at Supersonic Speeds" and presented at 
the July 1952 meeting of the Institute of Aeronautical Sciences. 
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/ 1.2 1.4 1.6 1.8 20 

Mo 

Sketch (e) 


Thickness Distributions Having Minimal Drag 
in Supersonic Flight 


The expression for the perturbation velocity potential given in 
equation (13) provides a ready /way of expressing analytically certain 
criteria that have been given by R. T. Jones (ref. ll) for wings pos- 
sessing minimum drag characteristics. An essential feature of Jones’ 
results involves the use of combined flow fields which are obtained by 
superimposing the disturbance fields in forward and reversed motions. 

So long as the governing equations of flow are linear, it is possible 
to establish reciprocity relations between the induced fields of arbi- 
trarily situated sources and doublets in combined flow fields. Condi- 
tions for minimum drag under imposed restrictions are then expressed 
in terms of the pressure induced in the superimposed fields. For 
example, it is found that if the thickness distribution for a symmetrical 
nonlifting wing yields a specified volume, then drag is a minimum if the 
thickness is distributed in such a way that the pressure gradient in the 
combined field remains constant over the plan form of the wing. The 
application of this condition to a slender wing in supersonic flight is 
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simple, since from equation ( 13 ) the perturbation potential in the two 
directions of flow can be written explicitly. In the plane of the wing, 
the forward flow yields 


^ A ^ Jn | y-35.1 dy, - |2 to a*. 


1 nb ^ dq^ 

^b 

and the reverse flow yields 


<p r <w°> - lr/ b a If 1^1 s s/Aj ln Ar 1 3X1 


where 2 is the streamwise length of the wing and b is local semi- 
span. Since 


d^P-p 

dn 


dtp. 


r 
dn 


the perturbation potential cpg in the combined field is 

-U 0 d ni 


cp 2 (x,y,0) = 


2 n dx 


/ S'CxJ 


2|z-xJ 

2 n dx x 

P 


and it follows directly that if S* ( 0 ) = S'(2) =0, pressure coefficient 
in the combined field is 


V(xJ 


r_ -iri2i 

*J 0 X-X X 


dx n 


( 38 ) 


In the case of thickness distribution with given volume, Cp^ is a 

linear function in x and equation ( 38 ) is precisely the same integral 
equation that arises in the determination of thickness distribution with 
given volume for a slender body of revolution in supersonic flight 
(refs. 12, 13, and l4) . The same chordwise distribution of area there- 
fore exists for wings and bodies of revolution under the given conditions. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 28, 1952.' ‘ 
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